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Muon catalyzed fusion in deuterium has traditionally been studied 
in gaseous and liquid targets. The TRIUMF solid-hydrogen-layer tar- 
get system has been used to study the fusion reaction rates in the solid 
phase of D2 at a target temperature of 3 K. Products of two distinct 
branches of the reaction were observed; neutrons by a liquid organic 
scintillator, and protons by a silicon detector located inside the target 
system. The effective molecular formation rate from the upper hyper- 
fine state of /id and the hyperfine transition rate have been measured: 
A3 = 2.7l(J)stat.{'i'2)sy3t. fJ.s~^, and A3 1 = 34.2(8)stat.(l)sast. ^J.s~^. 



*Present Address: Institut de Physique Nucleaire, Universite Catholique de Louvain, 
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The molecular formation rate is consistent with other recent measure- 
ments, but not with the theory for isolated molecules. The discrepancy 
may be due to incomplete thermalization, an effect which was inves- 
tigated by Monte Carlo calculations. Information on branching ratio 
parameters for the s and p wave d + d nuclear interaction has been 
extracted. 



I. INTRODUCTION 

For forty years it has been known that the introduction of a negative muon into a 
mixture of the three isotopes of hydrogen (protium, deuterium, and tritium) can lead 
to fusion reactions between the hydrogen nuclei, a process called muon catalyzed 
fusion or ^CF. The intricately connected molecular, atomic, and nuclear processes 
that occur have taken many years to identify and offer a richness of physics in both 
theoretical and experimental domains |H-H|- 

The TRIUMF solid-hydrogen-layer target system [|| was used to measure /zCF in 
solid deuterium for the first time. The preliminary results Mm from our experiment 
at low temperature were in sharp disagreement with theory and initiated further 
investigations of the reactions in solid D2 P-pl]| . This work presents our final results 
based on larger statistics and the simultaneous observation of two distinct branches 
of the fusion reaction. 

The /xCF reaction in pure deuterium has been well investigated for gaseous and 
liquid targets — see Refs. ||l^,|l3l and references therein — where the assumption that 
the D2 is not interacting with its neighbors has been used in the theoretical analy- 
sis. That approximation has led to generally good agreement between theory and 
experiment when the deuterium is a fluid. The data presented here from /iCF in 
solid deuterium maintained at 3 K indicate that it fails. 



II. THEORY 

A. Physical Processes 

A negative muon introduced into pure deuterium slows and captures on a 
deuteron to form an excited /irf in a characteristic time of picoseconds in high 
density targets such as solids or liquids |l4| - [l(| . The subsequent deexcitation of the 
muonic atom occurs via Stark, Auger, scattering, radiative, and transfer processes, 
which occur on the 100 ps time scale [|^. The process {i^d)ni +d ^i- {^id)nf +d with 
Uf < Ui, often called Coulomb deexcitation, provides a source of acceleration for the 
muonic atom since the final state contains only the massive particles which share 
the energy released in the Ui — > Uf transition. The subsequent energy distribution 
of the fid contains both thermal and epithermal (nonthermalized) components of 
muonic atoms. The relative pop ulations reflect the temperature and density of the 



muonic atoms' environment |18-2 



Once the atomic ground state has been reached, the reactions of the /id are 
dominated by the scattering from D2 molecules, including inelastic processes such 
as the realignment of the fid hyperfine state, or the formation of the molecular 
bound state dfid. 

A muonic deuterium atom may form a molecule or molecular ion via two mech- 
anisms. The first is the capture of the muonic atom on one nucleus of a deuterium 
molecule to form a molecular bound state with the subsequent release of the bound 
state energy by an Auger process: 

fid+ [ddee] -^ [{dfid)de] + CAugcr- (1) 

A second process, a resonance mechanism first proposed by Vesman Q, is a 
curiously fortuitous effect that depends on the energy levels of all the involved 
bound states. The existence of a state of the (dfid) bound by ^2 eV allows the 
excess energy of the system, the energy liberated in molecular formation plus the 
incident fid energy, to be absorbed in the excitation of rotational and vibrational 



states, {lyfKf), of the full six body final state (strongly bound states liberate energy 
above the dissociation threshold of the D2). A resonance in the formation cross 
section occurs when the incident energy of the iidp in hyperfine state F is such that 
the total initial energy of the system matches the energy of some set of final state 
excitation parameters, 

^dp + [ddee]^.^Ki "> [(df^d) j^ dee]^^j^^ . (2) 

Here, J and v refer to the orbital angular momentum and vibrational quantum num- 
ber respectively, while the ^'ij,and Kij refer to the vibrational quantum number 
and rotational quantum number of the six-body system. The resonant formation 
process has been well studied, resulting in calculations of the energy dependence of 
the rate H^. 

Calculation of an effective formation rate which can be compared to experiment 
must assume an energy distribution of the fid atoms in the D2 target. The dis- 
tribution, which is determined by the initial energy at which the fj,d was formed 
and the subsequent energy loss processes, is convoluted with the reaction rates at 
an appropriate temperature to produce an overall temperature dependent rate |I^ . 
The agreement between theory, which predicts a rapid thermalization of the fid and 
hence uses a fid population in thermal equilibrium with the surrounding deuterium, 
and experiments which use liquid and gas targets, is good (see Fig. m). The results 
of present measurements in solids are in disagreement. 

The small size characteristic of the muonic molecular bound states is such that 
nuclear fusion can occur. For deuterium, there are two fusion channels: 

dfid — > fi + t + p + 4.0 MeV 

fi + ^Iie + n + 3.3 MeV , . 

fi^Ee + n + 3. 3 MeY ^' 



Once fusion occurs, the muon is generally set free and can begin the process again 
(cycling); however, the Z — 2 helium nucleus can sometimes capture the muon into 
atomic orbital bound states stopping further reactions, a process called sticking. 
The number of fusion cycles in which the muon participates is limited by the muon 
lifetime and by the sticking probability. The kinetics for the cycle shown in Fig. ra 
are given in the Appendix. 

The energy difference between the two hyperfine states of the fid is 48.5 meV; the 
energy is readily liberated in processes which depopulate the upper spin state. There 
is a disagreement between theory and experiment for the hyperfine transition rate 
As 1 , a situation which is well discussed in Scrinzi et al. |12| . As indicated in Fig. g, 
theory models the shape of the temperature dependence, but the calculated rate 
is too high. Given the agreement between the calculated and measured formation 
rates (Fig. M and the accuracy in the calculations for the scattering reactions (which 
are expected to be correct to 10%), the discrepancy is not understood. 



B. Monte Carlo Studies 

The two-node approximation (Appendix) assumes that the rates used to repre- 
sent the various processes are time independent, or, equivalently, the result after 
convoluting the fid energy distribution with the energy dependence of the rates re- 
mains constant in time. Since the calculations of the muonic processes give energy 
dependent rates, it is natural to use the Monte Carlo method to study the time 
evolution of the muonic atom energy distribution and thus the interplay between 
all the energy dependent rates. 



Until recently, only scattering cross sections for niuonic atoms on free hydrogen 
isotopes were available. Monte Carlo simulations predicted that the time required 
to thermalize a muonic atom after formation was very short (20 ± 20 ns in D2 gas at 
5% of liquid density, shorter times at higher density). After rapid thermalization, 
the energy distribution of the atoms was constant in time making the formation 
rates time independent and thus validating the two-node approximation. 

At low energies, fid elastic scattering on a solid lattice does not reduce the /id lab 
kinetic energy due to the large mass of the recoiling crystal. It was initially suggested 
by M. Leon that mechanisms for energy loss might be curtailed for muonic atoms 
in solid hydrogen isotopes. First quantitative calculations of the energy dependent 
scattering cross sections in solids have been completed Q. In Fig. ^ the scattering 
cross section for F = | /id atoms on a solid ortho-deuterium target at 3 K and 
zero pressure is shown. The total cross section includes elastic scattering on the fee 
lattice of frozen D2 molecules and inelastic scattering processes — phonon transitions 
and rovibrational excitations of the D2 molecules. The ^d deceleration rate is 
determined by the inelastic processes because of the large recoil mass in elastic 
scattering processes. The rovibrational excitation threshold is about 10 meV, and 
below that the fj.d atoms lose energy in purely phonon scattering, an amplitude 
which also falls rapidly with decreasing fid energy. In Fig. the inelastic cross 
section for fid in solid D2 is compared to the doubled nuclear fid + d cross section 
which represents the first-order approximation to the free molecular D2 cross section 

US- 

Based on the new rate calculations, a Monte Carlo simulation predicting the time 
spectra of fusion products from an infinite solid D2 lattice has been performed. The 
resulting spectra were fit with the same kinetics function as the real data to extract 
effective rates. The analysis of the Monte Carlo spectra will be presented after 
discussing fits to the data. 

The suppression of the inelastic scattering probability in solid deuterium implies 
a slow thermalization rate for any epithermal fid atoms, which in turn means that 
the energy distribution of the muonic atoms can not be assumed to be constant in 
time. The implications of these effects on the use of the two-node approximation 
are given in the discussion section. 



III. EXPERIMENTAL EQUIPMENT 

A description of the cryogenic target system used in this measurement has been 
presented elsewhere P,pl|. In brief, a continuous-flow liquid helium cryostat was 
used to cool two 51 /xm thick gold foils to 3 K (Fig. ^. A gas spraying system 
directed the deuterium toward the adjacent surfaces of the foils where it condensed 
and froze, creating a solid deuterium film maintained in a high vacuum. A palladium 
filter, which operated above 600 K, was used to purify the D2 immediately prior to 
freezing so the ratio of even to odd molecular spins was expected to be statistical 
(nD2) |32|] . Muons entered the deuterium after passing through one thin gold foil, 
which acted as a degrader, and stopped either in the first or second solid hydrogen 
layer, or in the gold itself. Muon catalyzed fusion of dfid molecules then occurred. 
The higher energy fusion products, neutrons at 2.45 MeV and protons at 3.02 MeV, 
were detected. 

Measurements with solid layers of D2 and ^H2 were made. Two different 
layer thicknesses of solid deuterium, 10.54(53) mg • cm~^ (herein called thick) and 
2.11(11) mg • cm~^ (thin), ^ 40 cm^ in area, were frozen onto the foils pSJ. The 
thick D2 layer received an integrated muon flux of 220 x 10^ /i with 53(3)%^opped 
in the solid layer, while the thin layer received 164 x 10^ /i with a stopping fraction 



of 14(2)%. The study of the background in ^H2 (where there is no strong fusion 
process) for both thicknesses above and for bare target conditions was done with 
130 X 10^ /i. Table | summarizes the data. 



A. The Trigger 

The timing diagram for the trigger electronics is shown in Fig. ^. When the Tl 
scintillator (see Fig. ||) signaled the entrance of a muon, two gates were opened: B 
to indicate that a muon was in the target, and EVG, the event gate. If any of the 
detectors recorded an event at a time when EVG was true, such as the pictured TPn, 
the corresponding TRGn gate was opened. At the end of the event gate duration, 
the logical OR of all possible TRGn gates was placed in coincidence with the end of 
the event gate, EEVG, to determine whether the event was read (an event trigger) 
or cleared. Events for which a second muon arrived in the target within ±10 fis 
were discarded either by inhibiting the trigger or during analysis based on an event 
flag. An event trigger immediately initiated a hardware inhibit to prevent loss of 
the data already measured. A computer-driven signal maintained the inhibit while 
the data were read and while the CAMAC modules were cleared. An event clear 
would inhibit the trigger system and then remove any values accidently recorded 
by the CAMAC system. 



B. Neutron Detectors 

The arrangement of detectors is shown in Fig. g. Each neutron detector was a 
12.5 cm diameter by 10 cm deep cylinder of NE-213 liquid scintillator viewed by a 
photomultiplier. Hardware-based pulse shape discrimination was used to differen- 
tiate between photon and neutron events. Uncharged particles were selected by an 
anticoincidence with both members of a pair of charged-particle scintillators [EiNl 
and EiN2 for (i = 1, 2) in Fig. pj. Since photon events were very common, prescaling 
was done to keep the photon triggers to an acceptable level. Typical neutron rates 
during the measurements were less than 150 s~^. 

An effective method for background suppression in /iCF processes is the require- 
ment that the muon survive the event, a condition enforced by the detection of 
the muon decay electron. This delayed electron, or deLe, condition was satisfied if 
the muon decay electron was detected between 0.2 and 5.05 /iS after the time of a 
candidate fusion event. It strongly enhanced the signal-to-background ratio of the 
data. Electrons from muon decay were detected by the coincidence of either pair of 
charge- veto scintillators. 

The time spectrum of the decay electrons was used to determine the proportion 
of muons stopped in hydrogen. Decay electrons were generated either with a free 
muon lifetime (2.2 /xs), or the lifetime of muons in Au (79 ns). The intensity of 
each lifetime component was proportional to the number of muons stopped in the 
respective material. 



C. Silicon Charged-Particle Detectors 

A fully depleted, ion implanted silicon detector of 600 mm^ active area was 

mounted inside the vacuum system to detect the charged fusion products. Its 

150 /im active thickness was sufficient to measure protons of up to 4 MeV, and 



was well suited for detecting the fusion protons at (or below) 3.02 MeV. The detec- 
tor was mounted on the cryostat cold shield, in front of the neutron detector Nl, 
and viewed the solid deuterium layers directly (see Fig. 0). At the nominal operat- 
ing voltage of 30 V the best energy resolution achieved was 170 keV FWHM due to 
mounting constraints and the resulting long cable to the preamplifier. The 3 MeV 
protons produced by dd fusion were easily seen but the tritons and '^He were too low 
in energy and hence obscured in background or stopped in the solid target before 
reaching the detector. The rate from silicon detector events was about 25 s~^. 



IV. ANALYSIS 

A. Neutron Spectra 

Three main event types occur in the neutron spectrum: fusion neutrons, cap- 
ture neutrons, and neutrons from the ambient background (zone neutrons). Fusion 
neutrons, which constitute the signal, are generated with a time distribution Uf, 
according to the kinetics of the muon catalyzed fusion reactions in the target and 
are thus time correlated with the muon arrival and will satisfy the deLe condition 
with the efficiency Cg with which the electron detectors measure that a muon has 
decayed in the delayed time interval 0.2-5.05 /iS. This efficiency is comprised from 
a solid angle contribution and the time window efficiency: 

ee = r^eAo / dre-^^''. (4) 

Jo. 2 fj,s 

Neutrons from nuclear muon capture, with time distribution ric, are also time 
correlated with the arrival of the muon according to the capture and decay lifetimes, 
but are obviously not followed by a muon decay electron. Only those capture 
neutrons that are followed by some accidental event in the electron detectors will 
be admitted to the deLe spectrum (with efficiency ea). 

The neutron events from the ambient flux contribute a constant background, 
Uz- These events pass the deLe cut depending on the time t of the neutron with 
respect to the muon arrival to. The deLe condition requires an event in the electron 
detectors between 0.2 /is to 5.05 /is after the neutron. For a neutron that occurs 
more than 5.05 /iS before the muon arrival, only an accidental deLe condition with 
efficiency ea can pass the event into the deLe spectrum. As the time of the ambient 
background neutron event moves closer to the subsequent arrival of the muon, more 
of the cut interval overlaps with the time when a muon is known to be in the 
target and hence the efficiency for accepting the neutron grows. For a neutron 
event occurring 0.2 fis before the muon arrival, the time selection interval overlaps 
exactly with the muon arrival, hence the muon decay selects the background event 
with an efhciency slightly larger than the decay electron efficiency Ce- After the 
muon arrival, accidental events are accepted with the electron detection solid angle 
efficiency weighted by the exponential probability that a muon had survived to 
the time of neutron detection. The flat background, Hz, thus acquires a muon 
lifetime dependence in the delayed spectrum. This relationship is represented by 
the integral term in Eq. m) . The efficiency e^ appears in all terms since an accidental 
deLe condition permits neutrons into the deLe spectrum regardless of the neutron 
source. 

The raw (without the deLe condition) and deLe spectra are given by: 



V 



(5) 



and 

/ nt+5.a5f,s \ (6) 

V Jt+0.2 ps / 

where 0(r — tg) is a step function which models the discrete arrival time of the 
muon. 

Combining (||) and (||) gives: 

ndel.e = Ee"/ + f^aTlr 

r.t+5.05 fj.s \ (7) 



\ Jt+0.2 fiS J 



Note that knowledge of the unknown capture time spectrum, ric, is not required. 
The largest contribution to the background is from the accidental deLe coincidences, 
and is proportional to the raw distribution n^, a spectrum well measured during 
the experiment and having small uncertainties. The events in the deLe spectrum 
occurring in the short time interval measured before the muon arrival determine 
the value of n^ . 

In Eq. (|7|) there are two input spectra, each with an uncertainty. Replacing the 
integral by the symbol I{t), we have 

ndel_e ± Sndel_e = Ce'T-/ + Ca (?^r ± SUr) + I{t)n^. (8) 

The efficiency factor eg will eventually be absorbed in a fitted normalization con- 
stant, so it does not need to be determined independently. 

The weights for the x^ fit depend on the input spectra uncertainties; 

r 2 2i"^ 

weight = {5ndel_e) + {eaSur) . (9) 

The weights contain the variable Ca, so the fitting routine was iterated, each time 
re-evaluating the weights, until a stable value for 6^ was found. Figure M shows the 
data and the fitted curve. 



B. Proton Spectra 

The proton spectra exhibited many features, and careful study of five different 
types of target was necessary to extract meaningful parameters. The different target 
conditions were: bare target, thick and thin solid protium layers upstream and 
downstream, and thick and thin solid deuterium layers upstream and downstream 
(see Table i. 



1. Energy Spectra 

The fusion protons are produced at 3.02 MeV, but lose energy in the target 
resulting in a broad distribution: Fig. shows the spectrum taken with the thick 
deuterium target. The background comprises about 20% of the total counts. 

The background in the detector comes mainly from charged particle emission 
following nuclear muon capture. Protons are the most commonly emitted particle 
with a spectrum which is suppressed for energies below the nuclear Coulomb barrier. 
The emission probability ranges from 15% for muon capture on silicon — with a 



5 MeV Coulomb barrier — down to about 1% following muon capture on gold where 
the Coulomb barrier is about 16 MeV Q]. The relative intensity of the background 
depends on the thickness of the hydrogen layers and the incident muon momentum 
which together determine the ratio of muons stopped in the gold target support foil 
or detector silicon to those stopped in the hydrogen target itself. Thick hydrogen 
targets exhibited less background. 

The background in the silicon detector was removed to the e^ level by the deLe 
requirement, evidence that the background was due to a process which consumed the 
muon. The ea measured by the ratio of the number of counts above the 3 MeV fusion 
limit (hence nonfusion related), to those from the same spectrum after applying the 
deLe cut gave a value of 0.11(4)%, which is in agreement with the value found for 
Ca from the neutron detector: 0.14(4)%. 



2. Time Spectra 

Time spectra for each of the five target types were collected for 23 different 
energy windows, each 142 keV wide (the intrinsic energy resolution was 170 keV), 
beginning at 0.9 MeV. They were examined for number, intensity, and value of 
visible lifetimes, and time-zero value. For data where fusion events contributed to 
the spectrum the characteristic two-lifetime form of the kinetics in pure deuterium 
allowed an accurate measure of to for that energy window, as well as the intrinsic 
time resolution of the detector at that energy. 

The silicon detector exhibited time walk. Thick deuterium data (see Table g) 
were chosen to measure the energy dependent correction since the fusion protons 
provided a well defined time zero, and the energy loss in the layer distributed the 
signal over the energy range 1-3 MeV. Walk gave a 30 ns correction between the 
3 MeV and the 1 MeV time spectra, where the difference in proton flight times at 
the different energies was no more than 4 ns. 

Two standard energy regions were defined from which the time spectra were 
selected for further analysis. The lower energy signal region, called I, was limited 
to 2.3-2.7 MeV; below 2.3 MeV the time resolution was poor enough that it began 
to obscure the structure of the fast lifetime component in the fusion signal. The 
limit at 2.7 MeV was selected to favor fusions occurring away from the surface 
of the layer, where the possibility of ^d escape from the layer would not affect the 
kinetics representation. The higher energy background region, called h, was selected 
between 3.3-4.0 MeV, where there were no fusion events. Walk-corrected data were 
then selected for the h and I regions. 

Figure shows the time spectra for the h and I energy cut regions for the thick 
deuterium, bare target, and thick protium runs. Two lifetimes are the dominating 
feature in the thick deuterium data, (A), as expected from the two-node approxi- 
mation. The fusion spectrum has a rapid onset and rise of the leading edge, a shape 
given by the time resolution of the detector and not by the physical processes of 
muon stopping, atomic capture, and thermalization. 

The spectra without fusion signals, the h and I regions for both gold and thick 
protium and the h region for thick deuterium, seem delayed with respect to the 
fusion signal. It is possible to extract two lifetimes from the background spectra 
as well — lifetimes consistent with muons in gold (80 ns) and muons in intermediate 
mass nuclei like silicon or aluminum (~ 700 ns). The background does not pass 
the deLe requirement except by Eq coincidence, however the fusion statistics in the 
deLe spectrum were insufficient for the kinetics analysis. 

The background in the fusion spectrum was parametrized and removed phenom- 
enologically. Studies of the h and I time spectra for targets that did not produce 



fusion signals provided the method to predict the background in the / region given 
the spectrum in the h region. Scahng from the h region to the / region was time 
dependent, i.e., 

l{t) = f{t) X h{t). (10) 

This is consistent with a background of charged particles emitted after muon capture 
on differing nuclei. Each capture spectrum contributes to the overall time spectrum 
with an intensity dependent on the chosen energy range. Since the energy spectra 
of the emissions are nuclei dependent, the overall time spectrum cannot be expected 
to be independent of energy. The function f{t) ±5f, parametrized in four variables 
Ai via 

f{t) = [Ai + Aae-^^f*--^*) (l - e-^-^(*-^^))] , (11) 

was found using Eq. (HO) to fit the h and I time spectra from both thin and thick 
protium layers; 5f was found from the covariance matrix of the Ai. 

Since the deLe condition was not imposed on the fit data, the I spectra containing 
fusion signals were simply composed of the fusion protons and the background. The 
fit expression, with uncertainties, was 

Pi±5pi=pf + {f±5f){ph±5ph) (12) 

where pi was the measured time spectrum in the I energy window, ph was from the 
h energy window, pf was the kinetics function, after convolution with the detector 
resolution, for proton production by /^CF, and / was the background transformation 
function as above. 

The weight for each fitted point was constructed in the same manner as in the 
neutron case, yielding: 

weight = \{5pif + ifSphf + [phSff] "' . (13) 



C. Fit Methods 

The data were fit with a function made by convoluting the theoretically expected 
time distributions, Eq. (119) of the Appendix, with a Gaussian detector time reso- 
lution and a step function modeling the muon arrival. The time spectra data were 
rebinned with an adaptive step size which preserved the fine binning where sensitiv- 
ity to fast rates was required, but used larger steps at long times where small scale 
sensitivity was not needed. The larger bins at long times eliminated the problem of 
fitting data containing many zeros, and reduced the total number of points in the 
fit, which in turn significantly reduced the amount of computer time required per 
fit. 

Excluding background parametrization, six independent parameters were re- 
quired to completely determine the shape of a spectrum: two lifetimes, two am- 
plitudes, the detector timing resolution a, and the muon arrival time to. There are 
ten kinetic parameters which determine the values of the lifetimes and amplitudes, 
and so it is impossible to measure more than four of them at any one time. Since 
an absolute measurement was not done, one of the parameters was used for nor- 
malization, so only three of the ten kinetics parameters could be extracted from a 
single time spectrum. By fitting both neutron and proton data simultaneously, four 
kinetics parameters could be extracted (equivalent to an absolute measurement for 
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a single fusion product spectrum). The parametrization of the background in the 
data introduced other parameters which were independent of the kinetics. The fit 
to the data was made by x^ minimization. 

The fits to the spectra were sensitive to the values of As i, Ai, <j)zXz (the loss 

rate to target contaminants heavier than hydrogen, predominantly nitrogen), Ai, 
Ps, Ps, and Pp (see Fig. and the Appendix for the definition of the parameters). 
The standard values of the parameters passed to the fitting routine are given in Ta- 
ble O, along with the uncertainty on the value used when examining the systematic 
uncertainties. 



D. Fit Results 

For both neutron and proton data, fits were made to the individual time spectra. 
Neutron data were restricted to the thick deuterium targets only, since insufficient 
statistics existed for the thin deuterium layers. For the proton data, both thick and 
thin deuterium layers yielded spectra which could be fit. 

In the fits to individual spectra, the two kinetics parameters As i and Aa, as 
well as the loss rate to high-Z components, (pzXz, were measured. When the fits 
to the combined spectra were done, the sensitivity to four parameters was used to 
measure As i , As , (jiz^z, and values for Ai , Pg, Ps, and Pp, each taken in turn with 
the remaining three of the four fixed at the standard values (see Table |l|). 



1 . Uncertainties 

The systematic uncertainties fall into two categories; the effects on the measured 
values resulting from uncertainties in input parameters, and the effects due to vari- 
ations in the fit interval, cuts on the input spectrum, and, in the case of the proton 
spectrum, the background scaling function. The uncertainties due to input param- 
eters in the fit are relatively easy to identify and to understand by changing the 
parameter and refitting the data. The systematic effects due to the cuts on the spec- 
trum and the fitting interval were studied by fitting the same data with different 
cuts. The fit was considered satisfactory when the variations on the fit parame- 
ters with cut value were small with respect to the other systematic and statistical 
uncertainties. 

Within the stated uncertainty, the values found for the fitted parameters were not 
affected by variations in temperature T (used to make Ai 3 from the As 1 rate via 
detailed balance), density (j), muon decay rate Ao, the detection efficiency e (which 
was detector dependent), and Wg. This was checked by changing the input values 
of the above parameters by twice their respective uncertainties and confirming no 
variation in the fitted parameters. The uncertainty in the input values of Ai , Ps, 
Pp, and Ps did have significant systematic effects. 

The total systematic uncertainty on a fitted parameter was determined by the 
addition in quadrature of the uncertainties due to the possible variations of the 
fixed input parameters. The systematic uncertainty associated with the proton 
background scaling function was handled in a different way. The uncertainty in 
the background scaling function was explicitly taken into account in the fitting 
function, and so the uncert ainty due to the scaling is included in the statistical 



error evaluation (see § [V D 3 below) . Asymmetric uncertainties were determined by 
the change necessary to increase the x^ of the fit by 1.0, however, in no case were 
the uncertainties strongly asymmetric. 
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2. Neutron Data Fits 



Fits to the thick deuterium neutron data are summarized in Table |I| and plotted 
in Fig. n. The fitted value for the Eq efficiency agrees with values measured from 
neutrons above the fusion energy edge, measurements with the protium targets, 
and with measurements made with the silicon detector. The a value for the time 
resolution of the detector yields a FWHM of 3.5 ns, in accord with the flight time 
of 2.45 MeV neutrons traversing a 10 cm deep detector. 



3. Proton Data Fits 
Fits to the proton data taken from both the thick and thin deuterium layers are 



summarized in Table III ; the results of the fit for the thick target data are plotted 



in Fig. hOl The results for the thick deuteriur n gi ven in the table are in reasonable 



agreement with the neutron results in Table III. The signal-to-background ratio 
for the thick deuterium proton data was roughly 9:2, while the ratio for the thin 
deuterium data was three times worse at 9:6. The background scaling function 
was predicted using both the thick and thin protium layers and so, to check the 
systematic effects on the fits, the two deuterium spectra were each fitted with both 
of the scaling functions. The results of the fits showed that the thick deuterium 
data were not sensitive to the choice of background scaling, presumably due to the 
favorable signal-to-background ratio, while the thin deuterium fits fluctuated to 
the limits of the statistical accuracy. The data extracted from the thin deuterium 
targets were more sensitive to other effects such as /id escape from the layer and, 
from the value of (j)zXz, suffered from a higher contamination. 



4- Simultaneous Fits to Both n and p Spectra 

The simultaneous fit of both the neutron and proton data from the thick deu- 
terium target allowed one more kinetic parameter to be measured. The choice was 
limited to Ai, /3s, Pp, and Pg: the kinetics equations were not sufficiently sensitive 

to other parameters. All of the fits gave a x^ value of 325 for 309 degrees of freedom. 

With A 1 and any two of Ps, (3s, and (3p fixed, a value for the remaining parameter 

was found. The results of the fit, which yielded values for Pg, (3s, and (3p, in addition 



to All, Ai, and ((iz^z, are given in columns 2-4 of Table |IV . 

The values of As i. As, and (fiz^z do not change during the fits to Pg, (3s, and 
(3p, behavior which is expected from Eq. (Hq) which enters the kinetics as a distinct 
subformula. The large systematic uncertainty in the As value is due almost entirely 
to the uncertainty in the Ai input value. 

Fixing Pg, Ps, Pp at the standard values allowed a measurement of Ai, presented 



in the last column of Table IV. For this fit, systematic uncertainty in the parameters 



is dominated by the uncertainty in the Pp input value. 



E. Fits to the Monte Carlo Simulations 

A Monte Carlo (MC) simulation has been performed using the calcu- 
lated (Ref. 12^) energy-dependent rates for fj,d atom scattering on solid deuterium 
at 3 K, calculated rates of dfid molecule formation (resonant and nonresonant) on 
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a free D2 molecule and back-decay rates of the {dfj,d)d coraplex |25 2^. Since the 
solid target used was a statistical mixture of ortho- and para-deuterium (nD2), the 
/xd scattering on D2 molecules in both K = and K = 1 initial rotational states has 
been taken into account. The muon transfer from fj,d atoms to high-Z contamination 
has been described by a constant muon transfer rate. 

Recent experiments |£l|,Q and theory [|l^ show that the initial distribution of 
fid energy contains an admixture of nonthermalized (in the region of 1-10 eV) 
fid atoms. The initial fid energy in our MC has been described by a sum of two 
Maxwell distributions. One of them, corresponding to atoms thermalized at 3 K, 
has a mean energy of 0.4 meV. The hot fid atoms are assumed to have a mean energy 
of 1 eV. The details of the shape of the high-energy component are not important 
for this simulation because the deceleration rate at energies greater than about 
0.1 eV is very high. Therefore, the main input parameter to the Monte Carlo is the 
relative intensity of the energetic fid population. The second important parameter 
is a scaling factor for the energy-dependent nonresonant spin flip rate As i(£'), to 
account for the disagreement between theoretical results and experiment (see p^). 

A qualitative agreement with the two-node kinetics model has been reached by 
a proper choice of the two important parameters of the Monte Carlo. Fitted MC 
spectra yield values of As and A3 1 consistent with experiment if the scaling factor 
for the theoretical nonresonant spin flip rate is about 0.5 and the relative population 
of the hot fid's is near 0.75. However, this fit of the kinetics function to the early- 
time part of the neutron and proton MC spectra does not yield satisfactory values 
of x^- This result is probably due in part to the lack of a correct form of the 
energy-dependent dfid resonant molecular formation rate below about 10 meV, in 
solid deuterium at 3 K and zero pressure. 



Figure 11 shows the experimental neutron data in comparison to the output of 
the MC simulation. The spectrum at long times, when fid atoms disappear mainly 
via nonresonant dfid formation and muon transition to the Z> 1 contaminants, is 
well described by the MC simulations and good values of x^ ^re obtained (x^ = 
33.5/30 points). This is so since the rates for nonresonant dfid formation, in both 
the J = and J = 1 states, are not strongly energy dependent below 10 meV, 
meaning that the detailed shape of the fid energy distribution is not significant. The 
disagreement at times dominated by the quartet to doublet hyperfine transition and 
thermalization is clearly visible as the steeper slope in the MC spectrum. 



V. DISCUSSION AND CONCLUSIONS 

Results derived from our experiment show that the processes governing molecular 
formation in the solid state are not yet quantitatively understood, but qualitative 
consistency has been achieved with the theory of slow thermalization caused by 
solid state effects ||2^ . 



1. Kinetics Parameters 



The values presented in Table IV for the kinetics parameters (under the assump- 
tions inherent in the two-node approximation) represent the first examination of 
the muon catalyzed fusion process in solid D2. 

The qualitative success of the Monte Carlo simulation of fusion in solid D2 is a 
clear indication that slow thermalization effects are important. It is surprising that 
the two-node approximation works at all in this situation. Why it does work can 
be seen by examining the processes themselves. 
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For fids, slow thermalization at the lowest energies allows molecular formation 
via the lowest resonances. Even so, spin deexcitation processes and thermalization 
remove fid from the resonant energy region at least ten times faster than the effective 
formation rate. Under such conditions, sensitivity to the shape of the resonance is 
suppressed and use of a constant rate is adequate. 

Other explanations for the large A 3 rate again rely on the solid state of the deu- 
terium. Examination of Eq. dlj) shows that an increase in A/, the deexcitation rate 
leading to fusion for the resonantly formed [{dfjd)j^dee], will increase the effective 
formation rate by more successfully competing with the Tsf' deexcitation (back 
decay). If there exists a sufficiently strong coupling of the rovibrational modes of 
the resonant six-body complex to phonons in the solid, then there may be sufficient 
augmentation of the deexcitation rate to increase A/. The molecular formation 
process occurring in solids has not yet been rigorously examined. 

Since resonance formation involving fidi occurs at collision energies greater than 
50 meV, the rapid initial cooling of the fidi. permits only nonresonant molecular 
formation, thus assuming a time-independent nonresonant Ai rate is valid. Our 

measured value for Ai [— 0.052{8)stat.{^)sys. ^^s^^] is in agreement with nonreso- 
nant molecular formation rates measured in other experiments |36]. However, the 
assumption that Ai measured in gas is appropriate for 3 K solid contributes the 
dominant systematic uncertainty in our measurement of Aa . 

The Asi rate measures the disappearance rate for fid atoms capable of partici- 
pating in resonant formation via As and thus contains hyperfine deexcitation from 
direct scattering, a contribution from back decay, and the effective rate of energy 
loss for the fida which are in the energy window of the lowest resonance. Since our 
measured value for the effective hyperfine deexcitation is consistent with calcula- 
tions of the rate due to scattering alone, the calculated rates are too large, a result 
consistent with recent fid + HD experiments |36| . 

Emission of the fid from the finite solid target layer is a loss mechanism that is 
not explicitly accounted for in the analysis, however the process would mimic the 
(l)z^z loss mechanism. We believe that the contamination was due to nitrogen, based 
on the observation of muonic nitrogen x rays. The value found for 4>z^z, coupled 
with the value for the transfer to nitrogen, Xz ~ 10^ /^s^^, implies a nitrogen 
contamination of about 1.5 ppm, a value which was independently measured for 
the thick deuterium target via x-ray analysis. If the value were to represent only 
fid escaping from the layer, then roughly one third of all the fid atoms are emitted 
from the layer. Given the current understanding of the scattering cross sections, 
this process is unlikely to give such a strong effect. A Monte Carlo which includes 
the cross sections corrected for solid state effects as well as an accurate modeling of 
the finite geometry of the targets is required for the quantitative reinterpret at ion 
of the A3 and A3 1 rates in terms of thermalization. 

Demin et al. [ |ll[ have presented measured values for fiCF rates in solid deuterium 
over a range of temperatures, assuming a similar kinetics function. Their results 
are consistent with ours, thus the disagreement with theory indicates that further 
calculations for the molecular formation rate in solid hydrogen must be made. The 
next experiments will have to examine the thermalization process, as well as the 
possible para-ortho effects. One suggested increase in the molecular formation rate 
for thermalized fida occurs only for para-D2 {i.e., J = 1) plj . Using a solid ortho- 
D2 target would help to verify or discount that effect. 

In solid hydrogen the molecular formation rate is concentration dependent, pro- 
portional to Cd'/'No, however, the thermalization rate is proportional to 4>No, since 
the scattering and energy loss processes for both /irf+D2 and fj,d+^}i2 are similar for 
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low energy fid. If slow /jd thermalization is the process increasing the As molecular 

formation rate, an experiment in a solid mixture of D2 + ^H2 using several values 
of the deuterium concentration will not yield the same reduced formation rate A3 . 

Such an experiment in a D2 + ^H2 mixture would be more difficult to analyze, since 
the lip to fid transfer and subsequent creation of an epithermal fid would change the 
initial intensity of the epithermal fid atoms — the effect believed to be responsible 
for the elevated As rate. 



2. Branching Ratios 
Three values parametrize (3p^ the effective branching ratio: ps, Pp, and Pg 



[cf. Eq. (16)]. Table IV lists fitted values for each of these parameters determined 
using standard values for the other two. 

Our value Pp = 0.563(14)stat.(ll)sys. is close to the 0.59 predicted by theory 
(Hale, pql), and is consistent with a previous measurement by Balin et al. p9| . 
Our value, Ps — 0.487(15)sfat.(ll)sas., is consistent with 0.47, the prediction of 
Hale. 

Standard values for Ps and Pp were used in a fit to Pg since Ps and Pp, as discussed 
above, were consistent with both theory and experiment. Faifman p5[ predicts 
Ps = 0.560 at very low energy. Our result, Pg = 0.47(8)stat.(6)sys. is in agreement. It 
should also be noted that the Pg value is dependent on the energy at which molecular 
formation occurs, so for fid atoms in the process of thermalizing, a constant Ps does 
not accurately represent the true physics processes. 
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APPENDIX: REACTION KINETICS 

This appendix describes in more detail the kinetics equations derived from Fig. 0. 
All rates are normalized to liquid hydrogen density (4.25 x 10^^ atoms • cm^"^), and 
unit concentration, thus giving reduced rates for comparison with experiments under 
different conditions of concentration and density. The reduced rate of formation of 
dfid molecules from a fid atom in the hyperfine state F is Xp pq |: 

A^ = A„. + ^ Xps ~ ^^' ^ . (14) 

S ^/ + l^F" ^ SF" 

Molecular formation leading to fusion is composed of the nonresonant formation 
rate and the fraction of resonantly formed molecules which successfully fuse via A^ 
rather than resonantly scatter as characterized by Tsf' |2^ • Likewise, the reduced 
hyperfine transition rate from state F to F', Xff' , has contributions from regular 
scattering and from the decay of the resonantly formed [{dfid)dee] complex: 
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\FF' 



\PF' 



Ea 



FST 



TsF' 



^f + J2f" ^SF" 



(15) 



Finally, fip represents the fusion branching ratio for diid molecules formed from a 
fid in hyperfine state F. Experimentally, /?_f is given by 



(3f 



fid J. 



Do 



A« 



'He 



{fldF + B2^fJ. + p + t) + {fldF + B2 



>fj, 



mc 



The Pf measure will differ between the two hyperfine states due to the different 
distributions of bound state angular momenta selected by the hyperfine-dependent 
molecular formation processes p5| . The Pf parameters are convenient when writing 
the kinetics, but are composed of the more fundamental s and p wave fusion branch- 
ing ratios, together with Pg, the ratio of s to p wave bound states from nonresonant 
formation pM. Resonance formation always produces the p wave molecular state. 



f}F = ^ [FsPs + (1 - Ps) M + ^^^(3p 
Af ^f 



(16) 



The kinetics of the muon catalyzed reactions in pure deuterium are represented 
in Fig. y. Using effective rates, the time evolution of the two different hyperfine 
populations can be approximated as a differential equation of the form: 



d_fN3 



-A B 
C -D 



N3 
2 

Ni 



(17) 



with initial conditions rjF (normalized to a single muon) 



iV| (t = 0) 
7Vi(t = 0) 



Vi 



where 




The coefficients in the matrix can be read from Fig. ^ giving (here all explicitly 
positive): 



A = Ao 



02 Az 



{a. [l-,y.(l-e)(l-^3/3|)]+A..} 



i3 = 0[Ai773(l-6)(l-^3/3i)+Ai3 



C = 



k^il-^){l~^sPs) 



+ As 1 

2 2 



D = Xo 



'{h 



1 - r;i (1 - e) f 1 - u;,/3i 



+ A 



-I- 

2 2 J 



The formal solution of Eq. ( |l7| ) to obtain the time dependence of the hyperfine 
populations, NF(t), is tedious but straightforward. Once derived, the time structure 
of the fusion products, either protons or neutrons (generically denoted fc), can be 
found by summing over the populations, formation rates, and branching ratios: 



dk , Y^ X J.J .,, / Pf for k = n, 

—r — (P / aFAFJ^F with aF — { -, or? 
F=-.- 



(18) 
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The time distribution of product k is thus: 



dk ^ Lit ^ Lit 
—- = vl/ae 2 + ^le 2 

dt 2 2 

with rates which are the negatives of the the parameters: 
-1 



Ls = — 
5 2 



{A + D) + ^{A-DY+ABC 



(19) 



(20) 



and amphtudes: 



Li = 



*; 



(.4 + i^) - v/(^-^)'+4BC 



■^asAs 773IL3+Z)) 



773 IL3 +1) +7/1-B 



+ oiiXi 

2 2 



^'?l^-^^(^i+^)]} 



Li 



(21) 



(22) 



and, 



*i 



I as A3 

1 2 2 



+ ai Aj. 

2 2 



-r?3(L.+i?)-,?.B] 



(23) 



The total yisld of fusion particles from a single muon can be obtained by inte- 
grating Eq. ( [l9|) over times [0, 00) to obtain: 



Yk 




(24) 



The above representation of the kinetics is the exact solution to the two-node 
approximation used in Zmeskal et al. [nSl , and verifies the solution presented therein 
to within 1% accuracy. 



[2] 

[3] 
[4] 

[5] 
[6] 



[7] 



W. H. Breunlich, P. Kammel, J. S. Cohen, and M. Leon, Ann. Rev. Nucl. Part. Sci. 

39, 311 (1989). 

L. I. Ponomarev, Contemp. Phys. 31, 219 (1990). 

P. Froelich, Adv. Phys. 41, 405 (1992). 

J. S. Cohen, in Review of Fundamental Processes and Application of Atoms and Ions, 

edited by CD. Lin, World Scientific, Singapore (1993). 

P. E. Knowles et al, Nucl. Instrum. Methods A 368, 604 (1996). 

G. M. Marshall et al., in Muonic Atoms and Molecules, edited by L. A. Schaller and 

C. Petitjean (Birkhauser Verlag, CH-4010 Basel, 1993), pp. 251-260, [Proceedings of 

the Int. Workshop on Muonic Atoms and Molecules, Ascona]. 

P. Kammel, in Muonic Atoms and Molecules, edited by L. A. Schaller and C. Petitjean 

(Birkhauser Verlag, CH-4010, Basel, 1993), pp. 111-128, [Proceedings of the Centro 

Stefano Franscini, Ascona]. 



17 



[9; 
[lo; 

[11 
[12: 

[is; 

[14 

[is; 

[16; 

[ir 
[is; 

[19 
[20 

[21 
[22; 

[23; 

[24' 
[25^ 
[26' 
[27' 
[28' 
[29 
[30 
[31 
[32; 

[33 

[34' 
[35; 

[36 

[37' 
[38' 
[39 

[40 
[41 
[42' 
[43^ 
[44; 

[45' 
[46^ 



P. E. Knowles et ai, Hyp. Int. 101/102, 21 (1996). 

P. E. Knowles, Ph.D. thesis. University of Victoria, 1996, unpubhshed. 

P. Strasser, Ph.D. thesis. University of Tokyo, 1994, unpublished. 

D. L. Demin et ai, Hyp. Int. 101/102, 13 (1996). 
A. Scrinzi et al, Phys. Rev. A 47, 4691 (1993). 

J. Zmeskal et al, Phys. Rev. A 42, 1165 (1990). 

J. S. Cohen, Phys. Rev. A 27, 167 (1983). 

G. Y. Korenman and S. I. Rogovaya, Yad. Fiz. 22, 754 (1975), [Sov. J. Nucl. Phys., 

22, 389-393 (1976)]. 

M. Leon and H. A. Bethe, Phys. Rev. 127, 636 (1962). 

W. Czaphnski et ai, Phys. Rev. A 50, 518 (1994). 

E. C. Aschenauer et al, Phys. Rev. A 51, 1965 (1995). 
V. E. Markushin, Phys. Rev. A 50, 1137 (1994). 

J. B. Kraiman et al, Phys. Rev. Lett. 63, 1942 (1989). 

D. J. Abbott et ai, Phys. Rev. A 55, 214 (1997). 

E. A. Vesman, Pis'ma Zh. Eksp. Teor. Fiz. 5, 113 (1967), [JETP Lett., 5, 91-93 
(1967)]. 

L. I. Man'shikov et ai, Zh. Eksp. Teor. Fiz. 92, 1173 (1987), [Sov. Phys. JETP, 65(4), 

656-663 (1987)]. 

M. P. Faifman, L. I. Menshikov, and T. A. Strizh, Muon Catal. Fusion 4, 1 (1989). 

M. P. Faifman, Muon CataL Fusion 4, 341 (1989). 

A. Scrinzi, Muon Catal. Fusion 5/6, 179 (1990/91). 

M. P. Faifman et al, Hyp. Int. 101/102, 179 (1996). 

A. Adamczak, Hyp. Int. 101/102, 113 (1996). 

L. Bracci et al, Muon Catal. Fusion 5/6, 21 (1990/91). 

L. Bracci et al, Muon Catal. Fusion 4, 247 (1989). 

P. E. Knowles et ai, Hyp. Int. 82, 521 (1993). 

P. C. Souers, Hydrogen Properties for Fusion Energy (University of California Press, 

Berkeley, California, 1986). 

M. C. Fujiwara et ai, Hyp. Int. 101/102, 641 (1996). 

M. Lifshitz and P. Singer, Phys. Rev. C 22, 2135 (1980). 

F. Hartmann et al, Technical Report No. R-94-04, PSI, (unpublished), PSI Annual 
Report 1995/Annex I. 

C. Petitjean et al, Hyp. Int. 101/102, 1 (1996). 

L. I. Men'shikov and V.V. Filchenkov, Hyp. Int. 101/102, 207 (1996). 

G. M. Hale, Muon Catal. Fusion 5/6, 227 (1990/91). 

D. V. Balin et ai, Pis'ma Zh. Eksp. Teor. Fiz. 40, 318 (1984), [JETP Lett., 40, 
1112-1114 (1984)]. 

D. V. Balin et al, Phys. Lett. B 141, 173 (1984). 

C. Petitjean et ai, Muon Catal. Fusion 2, 37 (1988). 

D. V. Balin et al, Muon Catal. Fusion 5/6, 163 (1990/91). 
N. Nagele et al, Nucl. Phys. A 493, 397 (1989). 

V. P. Dzhelepov et al, Zh. Eksp. Teor. Fiz. 101, 1105 (1992), [Sov. Phys. JETP, 

74(4), 589-595 (1992)]. 

A. Adamczak, Muon Catal. Fusion 4, 31 (1989). 

P. Kammel et al, Phys. Lett. B 112, 319 (1982). 



TABLE I. Summary of the data. 



Hydrogen 


Name 


M 


% stopped in 


[mg-cm-2 -(foil)"^] 




(xlO^) 


hydrogen 


10.54(53) D2 


Thick Deuterium 


216.611 


53(3)% 


2.11(11) Da 


Thin Deuterium 


163.908 


14(2)% 





Bare Target 


18.504 


na 


1.06(6) 'Ha 


Thin Protium 


48.974 


14(2)% 
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5.27(27) ^Ha 



Thick Protium 



59.797 



53(3)% 





TABLE II. Standard values used for the fits. 


Parameter 


Value 


Source 


CJs 


0.122(3) 


Experiment; Balin et al. |4C(| 


d) 


1.4269(4) 


Experiment: R.H. Sherman in [if| 


Al3 

3 2 


O^s-i 


Detailed balance of A 3 i at 3 K 
- - 


Ai 

2 


0.044(5) iis"^ 


Review: Scrinzi et al. |12| 


P= 


0.560 


Theory: Faifmanj|25 


A 


0.47 


Theory: Hale, ||^ 


/3p 


0.580(5) 


Experiment: Balin et al. K] 



TABLE III. The kinetics values resulting from fits to the individual detector spectra. 



Parameter 


neutron 


pro1 


;on 




thick D2 


thick D2 


thin D2 


Ajx [ms-^] 


35.3(1.4) 


34.7(1.5) 


38.6(3.1) 


A3 [ms ^] 


3.12(14) 


2.77(12) 


2.97(23) 


^.A, [/is-^l 


0.288(15) 


0.340(12) 


0.589(24) 


fa 


0.14(4) % 


na 


na 


a [ns] 


1.48(11) 


9.52(28) 


9.02(36) 


xVdof 


158/150 


164/157 


162/157 


cl 


55% 


33% 


39% 



TABLE IV. The kinetics values resulting from a simultaneous fit to both neutron and 
proton spectra from the thick deuterium target, see the text for a complete explanation. 
The fixed values are taken from Table fn. 



Parameter 




Value ( Statistical) ( Systematic) 




Aji [^is-'] 




34.2(8)(1) 




34.0(8)(1) 


A3 [ms"^] 




2.71(7)(32) 




3.21(51)(16) 


4>z\z [ms~^] 




0.320(10)(1) 




0.320(10)(1) 


Al [ms-1] 


fixed 


fixed 


fixed 


0.052(8)(3) 


Ps 


0.47(8) (6) 


fixed 


fixed 


fixed 


A 


fixed 


0.487(15)(11) 


fixed 


fixed 


Pv 


fixed 


fixed 


0.563(14)(11) 


fixed 
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FIG. 1. Molecular formation rates found by the two-node approximation fit to data 
measured in solid, liquid and gas, and calculations of rates in liquid and gas. Statistical 
uncertainties are given, and the A3 values of Demin et al. KM have been normalized to 

Ai = 0.044 over their 5-30 K temperature range. 
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FIG. 2. Hyperfine transition rate Xsi. Theory calculations are shown for scattering 
only (dashed line) and for scattering plus back decay (solid line). 
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FIG. 3. The two-node approximation of the deuterium fusion kinetics. See the Ap- 
pendix for details and definitions. 
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FIG. 4. Calculated cross sections for scattering in D2. The dashed line shows the first 
order approximation to the cross section on free molecules. Note the logarithmic scale, 
the suppression of the inelastic scattering, and the Bragg cutoff for elastic scattering. See 
Pql for a description of the calculations. 
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FIG. 5. A top view of the detector positions with respect to the target. Not shown 
are the remaining two wire chambers, scintiUators, and Sodium Iodide crystal which con- 
stituted the imaging system, nor the Germanium detector, which, when in use, replaced 
neutron detector Nl. 
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FIG. 6. The timing diagram for the trigger electronics. Time is indicated horizontally, 
and NIM logic level on/off, indicated by the lower/higher lines. The relationship between 
the various signals is explained in the text. 
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FIG. 7. Neutron time data assiciated with Eq. m. Data points (A) are the deLe spec- 
trum. Points (B) (circle) show the ta-scaled contribution of the raw spectrum. Curve 
(C)(hne) is the representation of the other backgrounds, fit assuming a constant value 
before time zero, and a muon lifetime after time zero. The fitted contribution of the fusion 
neutrons (n/) is shown as the solid line passing among points (A). 
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FIG. 8. A plot of the silicon energy spectrum for solid deuterium. The 3 MeV proton 
peak is spread due to energy loss in the thick target. Note the strong suppression of the 
background in the deLe spectrum. The vertical lines show the h (dashed) and I (solid) 
energy cut regions. 
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FIG. 9. Silicon detector time spectra in the I (plus) and h (triangle) regions in A) thick 
solid deuterium, B) bare target, and C) thick protium target. For B) and C) the I spectrum 
is systematically lower than the h spectrum during early times. 
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FIG. 10. Curve (A) (plus) are the proton fusion data measured in the I energy region. 
Curve (B) (triangle) shows the background as determined from the h energy region and the 
scaling function f{t). The solid line is the fit and the dashed line shows the contribution 
of the fusion kinetics. 
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FIG. 11. The experimental neutron time spectrum (circles) in comparison with Monte 
Carlo predicted spectrum (line) scaled by a single normalizing factor. The inset shows also 
the uncertainty in the MC spectrum, illustrating the discrepancy at early times. 
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